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Abstract 

A particle system is a family of i.i.d. stochastic processes with values 
translated by Poisson points. We obtain conditions that ensure the sta- 
tionarity in time of the particle system in K d and in some cases provide 
a full characterisation of the stationarity property. It is shown how the 
characterisation problem relates to solutions of two-sided convolution 
equations. 

1 Introduction 

A Poisson process in the Euclidean space M. d is stationary if its intensity 
measure is proportional to the Lebesgue measure. More general Poisson 
processes can be denned on richer spaces, e.g. the space of functions or sets. 
While in these cases often there is no analogue of the Lebesgue measure, 
invariance properties of the process can be defined with respect to transfor- 
mations that account for the internal structure of the relevant phase space. 

One of most spectacular examples of this situation is due to Kabluchko 
[5], who considered the following situation. Let n be a Poisson point process 
in E and let > 1} be i.i.d. copies of a real-valued stochastic process 

£(t), t G R m . Define the family of functions x { + t G R m , for x t G EE, 
which (under appropriate integrability conditions on the intensity of the 
Poisson process) becomes a point process in the space of functions on M. m . 
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For any t € R m , N(t) = \x{ + £i(t) : i > 1} is the Poisson point process in 
R. Sometimes, the point process iV(i) formed by the values of the translated 
function is stationary in time even if £ is not stationary. It is important not 
to mix this concept with stationarity in R, where the points lie. 

Kabluchko [5] characterised the cases when a real-valued Gaussian pro- 
cess £ gives rise to a stationary point system N(t) called a stationary Gaus- 
sian system assuming that the intensity measure A of II satisfies J" R e~ £X A(dx) < 
oo for all e > 0. The stationary Gaussian systems are given by the following 
three classes. 

(i) A is an arbitrary measure on R and £ is a stationary Gaussian process. 

(ii) A is proportional to the Lebesgue measure on R and £(i) = W(t) + 
b(t) + c, where TV is a centred Gaussian process with stationary incre- 
ments, b(t) is an additive function, i.e. b(t + s) = b(t) + b(s) for all t 
and s, and c € R is a constant. 

(hi) The density of A is proportional to e~ Xx , x 6 R, with A ^ 0, and 
= W(t) — Xa 2 (t)/2 + c, where W is a centred Gaussian process 
with stationary increments and variance o~ 2 (t), and c £ R is a constant. 

The aim of this paper is to provide a partial generalisation of the above 
result for the case when £ takes values in a higher-dimensional Euclidean 
space, which is also mentioned in [5] as an interesting open problem. We 
show the relationship of the full characterisation problem to yet unsolved 
problems in real analysis. 

2 Multivariate particle systems 

Let {Ci,i > 1} be i.i.d. copies of a Revalued stochastic process £(t), t € 
R. All subsequent results can be easily generalised and remain valid for 
processes £ with argument t from a higher-dimensional Euclidean space. 

Furthermore, let II = {xj, i > 1} be a Poisson point process in R rf 
independent of the £1,^2, • • •• We call the process 

N(t) = {xi + &(t),i>i}, t€R, 

a particle system, so that a particle system is a stochastic process with 
values in the space of point configurations (or counting measures). Since 
the distribution of II is completely determined by its intensity measure A, 
we say that the particle system (A, £) is generated by measure A and process 
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By the finite-dimensional distributions of N we mean the distribution of 
the point process in R dn given by 

N(ti,...,t n ) = {(xi+&(t 1 ),...,Xi + Z i (t n )),i > 1}, h,...,tn€R. 

Denote by Pt x tn the finite-dimensional distributions of £, in particular 
Pt is the distribution of £(t). From now on we will always assume that the 
convolution A * Pt is a locally finite measure for all t G M. The following 
result is easy to obtain, e.g. using the probability generating functional of 
a Poisson process, see pp. 

Proposition 2.1. If A*Pt is a locally finite measure for all f £ R, then, for 
all t±, . . . ,t n € M., N(ti, . . . , t n ) is Poisson point process in M. dn with locally 
finite intensity measure 

A tl _ tn (A)= [ P tl _ tn (A-x)A(dx) (2.1) 

for all Borel A C M. dn , where A — x is A translated by (x, . . . , x) composed 
of n copies of x G 

The main question addressed throughout this paper is to characterise 
all pairs (A, £), such that the corresponding particle system iV is station- 
ary. Since the distribution of a Poisson point process is determined by its 
intensity measure, we immediately obtain the following result. 

Proposition 2.2. The particle system generated by A and £ is stationary 
if and only if 

M u ...,t n = Ati+s,...,t„+s (2.2) 

for all s,ti, . . . ,t n € M. 

3 Convolution equations 

The stationarity condition (|2.2p is in fact a system of convolution equations 
of the form 

P tl ,...,t n *A = P tl+s _ tn+s *A, (3.1) 

where A is the measure obtained by uplifting A to the diagonal in M. dn . In 
general notation, these equations are of the type 

o\ * fi = <72 * [i , (3-2) 
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where o\ and <7 2 are probability measures and \jl is an unknown locally finite 
measure. If 02 can be decomposed as <7 2 = o\ * o (or if o\ = cr 2 * a), then 
(|3.2p simplifies to 

(j, = /i * a (3-3) 

for another measure \i. This convolution equation was solved by Deny [2]. 
Namely, if the support of a is the whole R d , then all solutions of f|3.3[) are 
mixtures of exponential measures, i.e. 

H= [ t x Q(dX), (3.4) 
Je 

where t\ is the measure on ~R d with density e~( X ' x \ x G R d , and Q is a 
measure on the set E = E a with 

E a = {A g R d : / e< A ' a: V(dx) = 1} . (3.5) 

In particular, if £ is a real- valued Gaussian process with non-constant vari- 
ance cr 2 (t), t £ R, then there exist ti,i2 £ R such that c 2 (t2) > c 2 (ii) ; so 
that the first convolution equation Pt 1 * A = Pt 2 * A can be reduced to the 
Deny convolution equation fj3.3f) with a being the normal law of the vari- 
ance 2 {$2) — o 2 (t\). Hence A * Pf is necessarily a mixture of exponential 
measures, which is the crucial argument in the characterisation of stationary 
Gaussian systems in [5]. 

In the multivariate case the above argument does not apply any longer, 
since the difference of the covariance matrices of £(£2) an d (,(t\) may be 
neither positive nor negative definite. In the spirit of ()3.5|) . define 

E aia2 = {A G R d : f e^a^dx) = [ e^a 2 {dx)}. (3.6) 

JR d JR d 

While each measure fi given by (|3.4p with E = E aiCT2 satisfies the convolution 
equation (|3.2p . there exist solutions of (|3.2p not in the form (|3.4p . 

Example 3.1. Let &\ and o"2 be bivariate centred normal distributions with 
covariance matrices 

si= ( 1+ o Ci ?) and s2= G i+ci) • 

for some constants ci,C2 > 0. Let g : R — > R be a function such that 
f R g(x + y)e~ y l 2 dy is finite for all x G R. Then each measure [i g with the 
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density g(c x x\ + c 2 X2) satisfies (|3.2p . Indeed, substitution z = c 1 c 2 y 
yields that 

— / 5(— + — -— )e 2 ^ dy = — g(— + — )e H^J . 

ci 7k ci c 2 ci c 2 ,/ k c x c 2 c 2 

It remains to note that the both sides of this equality are up to the same 
constant the densities of the convolution of \i g and the centred normal dis- 
tributions with covariance matrices 

(0 0) and (0 ci) 

Note that E aia2 = {(Ai,A 2 ) : ci|Ai| = c 2 |A 2 |}. For instance, if g(x) = x 2 
then there exists no measure Q such that fj, g = f E e\Q(d\). 

Unfortunately there is no general result describing solutions of (|3.2p . 
The two-sided convolution equation can we written as /x * v = for a signed 
measure v with finite total variation. If v has a bounded support, then the 
density of /x solving this equation is called mean periodic function. Typical 
examples of mean periodic functions are exponential polynomials, i.e. sums 
of products of polynomials and exponential functions. While exponential 
polynomials are dense in the family of mean periodic functions on the line 
[3], this is unknown for higher-dimensional spaces. The situation with v 
having unbounded support (like in case of Gaussian measures) is even less 
explored. Thus, there are no analytical tools suitable to identify solutions 
of the two-sided convolution equation (|3.2p . 

4 Multivariate stationarity 

In this section we characterise the stationarity conditions for some (but still 
rather general) intensity measures A. 

4.1 Exponential measures. 

Consider candidates for the solutions of (|3.ip of the form A = for A G M. d . 
It is easy to see that necessarily A S Ep t p g (see (13. 6p ) for any t, s € R. The 
convolution e\ * Pt is locally finite if and only if 

Ee {x ' m < 00 for all t G R. (4.1) 
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Then the characteristic function with a complex argument in its first coor- 
dinate 

^i,...,t n («l -i\,...,Un) 

= Eexp{t(((«i - »A),^(ti)) + (u 2) £(i 2 )) + • • • + («n,C(tn)>)} 

exists for all «i, . . . , u n G M d , where i is the imaginary unit. 



Theorem 4.1. Assume that ^4-l\ ) holds. The particle system generated by 
e\ and £ is stationary if and only if 

<Pti,... ) t n («i - «A, ... ,u n ) = <pt 1+s ,...,* n + a (ui - ... ,Un) (4.2) 

/or all n > 1, s, t\, . . . , i n € K and tti, . . . , u n € M. d satisfying Y^i=i u i = 0- 

Proof. The proof follows the idea of [6j Prop. 6]. Let A be a bounded Borel 
set in R dn . Then 

K,..,tn( A ) = E / 1 («ti)+x,...,C(t n )+a:)e>i e-< A '^dx 

^.....tnC^-^e-^cfa, (4.3) 



where [i is a measure on M rfn given by 



^,...,* n (A) = E 



L (°.*(*a)-«ti),...,«(t»)-eOSi))eA e 



<A,f(*i)> 



(4.4) 



Since \it\,...,t n is supported by the subspace {(xi, . . . ,x n ) G R dn : xi = 0}, 
the decomposition A il j f n ( A) = f Hti,...,t n ( A ~ z)e~( x,z ^dz is unique, e.g. see 
Th. 15.3.3]. Finally, note that the Fourier transform of Ht x t„ is given 

by 



Vtl,...,t n (Ul, ...,U n )= ip tl ,...,tr, 



»A - y^Uj,u 2 , ...,u r , 

i=2 



□ 



A similar proof with the Laplace transform instead of the Fourier trans- 
form yields the following result. 

Proposition 4.2. Assume that the Laplace transform 

^ ll ...,t„(«l, • • • ,Un) = Eexp{(ni,£(ti)) + • • • + (u n ,£(t n ))} (4.5) 

exists for all u\, . . . ,u n £M. d such that Y27=i Ui = ^- Then the particle sys- 
tem (ca,0 is stationary if and only ifi>t x ,...,tjux, ...,u n )= ip tl+s> ... jtn+s (ux, ...,u r 
for alln> 1, s, ii , . . . , t n 6l and u\, . . . , u n € M. d satisfying 5^" =1 u\ = A. 
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For Gaussian processes we can give a more precise statement. Denote by 
E(ti,i2) the covariance matrix of and £(£2), in particular E(t,i) is the 
covariance matrix of £(t). It is important to note that, unlike the univari- 
ate case, S(ti,t2) may differ from E(i2,ii), namely E(t 2 ,ii) = E(ii,i2) . 
Furthermore, the covariance matrix (variogram) of £^2) ~~ is given by 



We say that multivariate Gaussian process £ has wide sense stationary in- 
crements if and only if T(ti,t2) depends only on the difference t\ — ti. In the 
univariate case, this property is equivalent to the fact that £(t + s) — 
t £ R, is stationary for each sGl, see Lemma 1], while in the multivari- 
ate case this is not so. 

Example 4.3. Let ^(t) = W{t) and let £ 2 {t) = W(t + h) for some fixed h, 
where W is the Wiener process. Then E£ 1 (ti)£ 2 (t2) is not necessarily equal 
to E£ 1 (t2)C 2 (^i)i so that E(ti,t2) is not necessarily symmetric. 

Theorem 4.4. TTie measure e\ together with a Gaussian process £ generate 
a stationary particle system if and only if 



where W is a centred Gaussian process with wide sense stationary increments 
and variance E(i,i), c S R rf is deterministic, and b : R — > M rf is o function 
orthogonal to A suc/i i/iai 



depends only on the difference ti —t\. 

Remark 4.5. If A = 0, condition (|4.7p implies that — 6(0) is an additive 
function, see (5J Lemma 2]. This is also the case if S(ii, £2) is symmetric for 
all ii and £2, e.g. in the univariate case where the orthogonality of b and A 
implies that b vanishes if A ^ 0. 

We use the following lemma, that is is easy to prove by a direct compu- 
tation. 

Lemma 4.6. Consider all Gaussian vectors in the Euclidean space W 1 whose 
Laplace transform ip{u) is given for all u from L + a, where L is a linear 
subspace ofW 1 and a £ R n . Then all these vectors share the same values of 
A T T,A, A T (m+Ea) and (rn, a) + |(o, Sa}, where m andY* are the mean and 



covariance matrix of the corresponding vector and A denotes any projection 
ofW 1 onto L. 



T{t 1 ,t 2 ) = E(t 2 , t 2 ) - E(ii, t 2 ) - E(t 2 , h) + E(ti, h) . 




(4.6) 



- b(h) + -(E(t 2 ,ti) - S(t!,t 2 ))A 



(4.7) 
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Proof of Theorem \4-4\ The sufficiency follows by the explicit writing of the 
Laplace transform of (£(ti), ■ ■ ■ , £(t n )). For the necessity, let n = 2 and apply 
Lemma B~6l and Proposition 14.21 with L = {(^1,^2) € M 2d : u\ + u% = 0} 
and a = (A, 0) € R 2d . Define block matrices 



A 



I -I 
-I I 



'E(ii,ii) S(t!,t 2 ) N 
,£(t 2 ,ti) E(t 2 ,t 2 ), 



where 7 is the d-dimensional unit matrix. Note that A defines a projection 
on L and E is the covariance of (f(*i)>£(*2))- Then all elements of ,4 T X;,4 
are proportional to T(ti,t2), meaning that W(t) = £(t) — E£(i), iel, has 
wide sense stationary increments. 

Define m(t) = E£(i). Calculating A J (m + So) with m = (m(ii), m(i2)) 
it is easy to see that 

m(t 2 ) - m(ti) + (E(ti,ti) - S(t 2 ,ti))A (4.8) 

is invariant after (ti,i 2 ) is replaced by (ii + s,t 2 + s). Denoting 

6(t) = m(t) + 2 E (t,t)A 

and using the fact that r(ii,i 2 ) = T(ti + s,i 2 + s), we arrive at (14. 7p . 
Furthermore, 

(m, a) + -(a, Ea) = (m(ti), A) + -(A, E(ti, *i)A) = (b{t), A) 

does not depend on £1, so that (6(i), A) is constant. Finally, set c = 6(0) and 
replace b(t) by b(t) - 6(0). □ 

4.2 Mixtures of exponential measures. 

Assume that A = J E t\Q(dX), where Q is a measure supported by E C 
so that A is locally finite. 

Theorem 4.7. Assume that (|4.ip ZioZds /or a// A /rom an open neighbour- 
hood U of E. The particle system generated by A = J t\Q(dX) and £ is 
stationary if and only if the system (t\,£) is stationary for all A € E. 

Proof We only need to prove the necessity. For w£l' i define 

E x = {A € E: {X,v) < 1}, 
E 2 = {\eE: (X,v) > 1}. 
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Let Aj = f E . t\Q(d\), i = 1,2. Without loss of generality assume that 
neither Ai nor A2 is the zero measure. Let A be a bounded Borel set. Since 
A satisfies (1231) . 

Ai * (P tl ,...,t n - P tl + s ,...,t n +s)(A) = A 2 * (P tl+s ,...,t n+ s - Pt^tJiA) . (4.9) 
Assume that (|4.9p is positive. Since Q{E{) > 0, 

Ai * (Pt u ...,t n - P tl +s,...,t n +s)(A + v)> * {P tl ,...,t n - P tl + S ,...,t n+S )(A) , 

A 2 * (Pt 1+S) ...,t n +s ~ Pt^tJiA + v)< e~ l k 2 * (P tl+S ,... )tn+S - P tu ...,t n ){A) . 

In view of (|4"Uj) . 

Ai * {Pt-i,...,t n - P tl + s ,...,t n +s)(A + v)>A 2 * (Pt 1+s ,...,t n +s ~ Pt u -,t n )(A + v) . 

Rearranging the terms yields that 

(Ai + A 2 ) * P tl ,...,t n (A + v)> (Ai + A 2 ) * P tl+St ... u+S (A + v) , 

which contradicts that A = Ai + A 2 satisfies (|2.2p . A similar argument 
excludes the negativity of ()4.9|) . and therefore Ai and A 2 satisfy (j2.2j) for all 
bounded Borel A. 

Consider any Ao £ E. By cutting E with hyperplanes, it is possible 
to construct a sequence of relatively compact sets Ef. C E, k > 1, such 
that Ek 4 {Ao}, the closure of E\ is a subset of U and A& = J E t\Q(d\) 
satisfies (|2.2p for all k. Since (|2.2p is scale invariant, it also holds for A/, = 
J Eh txQk(d\) with Q k (-) = Q(-)/Q(E k ). For all k, 

inf e~ {x ' x) < [ e~ {X ' x) Q k (d\) < sup e~ {X ' x) , (4.10) 
A e£fc JE k xeE k 

Since the both sides of (|4.10p converge to e~( A °'^ , A fc (A) ->• e,x (A) for all 
measurable A. 

It remains to show that the limiting measure satisfies (|2.2p . By (|4.3|) . 

l k *P tl ,...,t n (A)= f [ ^ tl ,...,t n (A-x)e-^Q k (d\)dx, (4.11) 

where f-ki,...,t n (A) is defined in (|4.4p . 

Since Qfc(^ fc ) = 1 and - x) < 1^ (x)Ee< A ' f (* 1 ^ , where A 1 is 

the set of x £R d such that (x, y) £ A for some y E M d ( n_1 ), 

/ Mt lr .. lt n(^ - xJe-^QfcCdA) < cl Al (x)e-^ , 
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where c is the supremum of Ee^'^* 1 " for A from the closure of E\. This 
supremum is finite, since Ee^'^*^ is analytic, hence continuous, in its do- 
main U. Since A\ is bounded, 

I cl Al {x)a' {Kx) < oo 

and the Lebesgue dominated convergence theorem yields 

lim A k *P tl ,... ttn (A)= I lim / Ht x ,... u {A - x) e -^Q k (d\) dx 

k^oo J R d fc->oo J Ek 

= I »t u ...,t n (A - x)e~^dx = e Ao * P tu ... ttn (A) , 

where the second equality follows by a simular argument as (|4.10p . □ 

Thus, if A is a mixture of exponential measures, then the mixing measure 
Q is supported by the set Dt lt t 2 mEp tl p t2 ■ 

Proposition 4.8. Let Ai, A2 € M. d with Ai 7^ A2. If the Gaussian systems 
(eAuC) an d ( e A 2 7^) are stationary, then the one- dimensional stochastic pro- 
cess (£ — E£, A2 — Ai) is stationary. 

Proof. Writing (|4.6p for Aj, i = 1,2, we arrive at 

W(t) - is(t, t)Aj + 61 (t) + ci = W(t) - ^E(t, t)A 2 + 62 (*) + c 2 • 
Since £(i 2 ,ii) — ^2) is a skew symmetric matrix, (|4.7p implies that 

(A 2 ,6i(t 2 ) - 6i(*i)> + (Ai, 6 2 (t 2 ) - &2(*i)) (4.12) 

is invariant after (ti, <2) is replaced by (ti + s, £3 + s). Denote shortly AA = 
Ai — A 2 . Rewriting (j4.12j) yields 

(AA,E(ti,ti)AA) - (AA,S(t 2 ,t 2 )AA) 

= (AA, S(ti + s, h + s) AA) - (AA, S(t 2 + s, t 2 + s) AA) . 

By [H Lemma 2], the function (AA, £(i, t)AA) is an additive function plus 
a constant. In view of the positive definiteness of T*(t,t), we conclude that 
(AA, £(t, t)AA) is constant for all t. The statement follows from the fact 
that a Gaussian process with stationary increments and constant variance 
is itself stationary. □ 
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The following result characterises stationary particle systems in case the 
two-sided Deny equation reduces to the one-sided one. 

Corollary 4.9. Assume that £ and A generate a stationary particle system, 
where £ is a Gaussian process such that P^ = Pt 2 * a for a Gaussian measure 
a and some t\^t%. If no linear combination of the components of £ — E£ 
is stationary, then A = ct\ for some c > and £ is given by (|4.6|) . 

Proof. The Deny theorem implies that A is a mixture of exponential mea- 
sures, so the result follows from Proposition 14.81 □ 

In particular, Corollary 14.91 applies if is a.s. deterministic for at least 
one t, for instance if £(0) = 0. Furthermore, it yields the result of [5] for 
non-stationary univariate process £. 

Example 4.10. Let ^(t) = £ 2 (t) = W(t) - at/2, where W is the standard 
Brownian motion and a € R. Then A = J K t( a +\-\)Q{d^) for a measure 
Q on R satisfying the integrability condition and £ = (£ 1 ,£ 2 ) generate a 
stationary particle system. 

4.3 Measures with exponential polynomial densities. 

Assume that A has the density 

p(x)e- {x ' x) = c a x a e { ~ x ' x) , 

\a\<k 

where p(x) is a non-negative polynomial of degree k. We use the multi-index 
notation, i.e. a = (a 1 , . . . , a d ), \a\ = a l + - ■ -+a d and x a = (x 1 )" 1 • • • {x d ) ad . 
Note that one can also consider solutions of convolution equations with not 
necessarily non-negative polynomials, which however do not admit a point 
processes interpretation. Nonetheless, even then we speak about stationary 
particle systems. 

Theorem 4.11. If the particle system (p(x)e~^' K ' x \^) for a polynomial p 
is stationary, then the particle system (q(x)e~^ X ' x \^) is stationary for each 
polynomial q obtained as a partial derivative of p. 
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Proof. For each n, bounded Borel set A in M. dn and x £ R d , 
A tl ,..., tn (^ + ^)= / Pt u ...,tM + x-z)p(z)e-^dz 

Pn,...,t n {A - u)p{u + x)e~^ u+x Uu 
= E^ e " <A ' X> / P tl ,...,t n (A-u)qp(u)e-^dn, 

/3>0 P ' jRd 

where qp is the partial derivative of p of order /3 and /3! = /3 1 ! • • • /3 d !. The 
stationarity of the particle system and the uniqueness of the polynomial 
imply that the coefficients of the polynomial do not change, and so the 
statement of theorem follows. □ 

Theorem 4.12. The process £ and A with density Y17=iPi( x ) e ~^ Xi ' X ^ 9 en ~ 
erate a stationary particle system if and only if £ and measure with density 
Pi(x)e~( Xi ' x } form stationary particle systems for all i = 1,... ,n. 

Proof. As in the proof of Theorem 14. 7} define E\ and E 2 , where now both 
these sets are finite. Since the exponential grows faster than polynomial, for 
each bounded Borel A and sufficiently large h, 

Ai * (P tl ,...,t n - P tl+s ,...,t n +s){A + hv) > e- x Ai * {P tl ,...,t n ~ Pt 1+s ,...,t n +s)(A) , 
A 2 * (P tl+s ,...,t n +s ~ P*,...,tn)(A + hv) < e~ l A 2 * (P tl+s ,...,t B+a - P tl ,...,t n )(A) , 

which eventually leads to a contradiction as in the proof of Theorem 14. 71 □ 

Theorem 4.13. Assume that Ee^'^*^ < oo for all t 6 R and all u from 
an open neighbourhood of X. Then the process £ and the measure with expo- 
nential polynomial density p(x)e~ ( - X ' x ^ generate a stationary particle system 
if and only if 

d n 
qi-Q-^ti,...,^-™ ~ ^2ui,u 2 , ■ ■ .,u n )\ x=x 

X i=2 

d - 
= Q(Q-) c Pti+8,...,t n +s(-' lx -^2ui,u 2 , ■ ■ .,u n )\ x=x , (4.13) 

X i=2 

for all partial derivatives q of p, all n > 1, s,t±, . . . ,t n € R and u\, . . . , u n € 
R d with Y%=i Ui = 0. 
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Proof. Denote shortly A£ = (£(t 2 ) - 
L~3lh for a bounded Borel A, 



>C(*n) - Similarly to 



A, 



E 



1^(0, A£)e< A ^)>p(* - e~< A '^ 



/3>0 



E 



R d 



LA-z 



0Jx£(.h)) 



x=X 



i{z)e~ {x > z) dz 



/3>0 17 K 



n (^-z)U =A(?/3 (z) e -^dz, 



where qp denotes the /3'th partial derivative of p. By Theorems 14. II and 14,111 
the value of tn(-A) is invariant for time shifts if and only if all the 
partial derivatives of \i are invariant. Taking the Fourier transform yields 
the claim. □ 

Now assume that £ is Gaussian. If £ and A with density p(x)e~^ X ' x ^ 
generate a stationary particle system, then £ and e\ also do, so that £ is 
described by Theorem [ 



Example 4.14. While in the univariate case Gaussian systems with a posi- 
tive exponential polynomial density do not exist unless the polynomial part 
is constant, the convolution equation can be satisfied with a signed measure 
A. For instance, the one-dimensional signed measures with density x 2k+l , 
k > 1, together with the standard Brownian motion form stationary particle 
systems. 



4.4 Exponential measures on subspaces 

Now assume that A is supported by a linear subspace EI of M d . Denote by 
the measure on H with density e~( X ' x ' , x & EL The corresponding Poisson 
point process is then a subset of EI. Without loss of generality, it is possible 
to assume that A € E and otherwise consider its orthogonal projection on H, 
which results in the same density. Denote by £ (i) the orthogonal projection 
of £(t) onto EI and let ^ = £ - £ H . 

Theorem 4.15. Assume that (]4.ip holds. The process £ and gener- 
ate a stationary particle system if and only if (|4.2p holds for all n > 1, 
s,ti, . . . ,t n 6l and u\,. . . ,u n € M. d such that Yli=l Ui ^ s orthogonal to E. 

Proof. The proof is similar to the proof of Theorem 14.11 with 



/i(A) = E[l ( £ (tl) _£H (tl)i .. .,£(t n )-^ (tl))eA e <A ' 5 (tl)> ] 
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□ 



The following theorem concerns the Gaussian case. Let m , m be the 
expectations of £ H , ^ _L . Furthermore, let E H (ii,i 2 ) (respectively E (ii,i2) 
and C(ii,i 2 )) be the covariance matrix of £ H (ti) and £ H (t 2 ) (respectively of 
^{h) and ^(t 2 ) and of £ m (h) and ^(t 2 ))- Finally r H (ii,i 2 ) denotes the 
variogram of £ H . 

Theorem 4.16. A Gaussian stochastic process £ and measure generate a 
stationary particle system if and only if £ satisfies the following conditions. 

(i) £ H has representation (14. 6p described in Theorem \4--4\ 

(ii) £ — is stationary. 

(Hi) C(t\ + s, ti + s) = C(ti, i 2 ) + -B(s) /or aZZ s, ii, i 2 £ Rj where B(s) is 
a matrix valued additive function. 

(iv) m L {t) = m-^O) - £(t) T A for all t € R. 

Proof. By applying a linear transformation, it is easy to reduce the situation 
to the case of A supported by the plane M spanned by the first k < d basis 



vectors in R d . Iff(t) = (^(t), . . . ,£ d (t)), then £ H = (^(i), . . . , £ k (t), 0, . . . , 0) 
and ^ = (0,...,0,£( fe+1 )(i),...,£ d (i)). By Theorem HJ22 consider the 



Laplace transforms with A — ^ U{ being zeroes in its first k coordinates. As 
in Theorem 14.41 consider the space L that contains (u\, u 2 ) with u\ + u 2 = 
for % = 1, . . . , k and a = (A, 0). Then 



is the projection on L, where Ik is the matrix with first k diagonal entries 
being one and otherwise zeroes. Then 




A 1 XA 




~ C22 



Cn — C21 

^11 

C21 — C11 




L/ 22 u 21 




(4.14) 



/m? + S=A - m? - S=A\ 



A T (m + Ea) 



(4.15) 



(m, a) + I (a, Ea) = (mf , A) + ~<A, E= A) , 



(4.16) 
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where £y = Y>(ti,tj), = T(ti,tj) and Cij = C(ti,tj). The invariance of 
r?j, the first row of (|4,15p and ()4.16j) imply the representation of £ H . The 
invariance of in (|4.14p yields the stationarity of — m 1 - . 

Define Bt 1; t 2 ( s ) = C(ti+s,t2+s)—C(ti,t2). The invariance of C(ti,t 2 ) — 
C(ti,t%) and C(t2,h) — C(t\,t\) with respect to translation of the time ar- 
guments yields that B tut2 (si + s 2 ) = B tltt2 (si) + B tut2 (s2) and that B tljt2 (s) 
does not depend on £1, £3- For every matrix valued additive function B the 
covariances of the form C{t\ + s,t2 + s) = C(t\,t2) + B(s) satisfy the in- 
variance properties. Finally (iv) is obtained by considering the second and 
fourth row of (|4. 15 j) . 

For the sufficiency note that the Laplace transform of the random vector 
(£(ii), • • • , at the point 




consists of a combination of similar elements as given by (|4.14p . (|4.15j) and 
(|4.16p . where and uf~ denote orthogonal projection of u- L on H and its 
complement. □ 

The following example shows that there exists stationary systems gen- 
erated by a process where £ is not stationary. 

Example 4.17. Let I = 1 x {0} C M 2 and consider £ = (£\£ 2 ) with 
^(t) = Zt — \\t 2 and £ 2 (i) = Z — Xt for the standard Gaussian variable 
Z and A € M. By Theorem 14.161 6 together with ^ form a stationary 
particle system. 



5 Multivariate Brown— Resnick processes 

Consider a special case of the particle system that appears if the Poisson 
process II lives on the diagonal line H = {x 1 = ■ ■ ■ = x d } in M. d . In this 
case, instead of the additive particle system it is convenient to consider the 
multiplicative particle system 

N e (t) = {y l e^ t \i > 1}, t€R, (5.1) 

where {yi : i > 1} = II e is a Poisson process on (0, 00) with intensity 
measure A e and independent of n > 1}, which are i.i.d. copies of a 
HL d - valued stochastic process satisfying 

Ee m < 00 foralUeK. (5.2) 
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Note that the exponential is applied coordinatewisely and the finiteness of 
expectation means that all its coordinates are finite. Then the intensity 
measure of N e (ti, . . . ,t n ) is locally finite and given by 

K,.,tM)= \ P{(e f H...,e^)ey-^4}A e ( ( i 2 /) (5.3) 

for all Borel A C R dn . 

Assume that IF has intensity measure A e (dy) = y~ 2 dy and define a 
process 77 with values in R rf by 

00 

rj(t) = \J yiffr®, teR, (5.4) 
i=i 

where the maximum is taken coordinatewisely. It is well known that the 
process 77 is max-stable with unit Frechet margins, see [3]. In order to deter- 
mine the finite-dimensional distributions of 77 note that the event {r](ti) < 
z\, . . . ,rj(t n ) < z n } (with coordinatewise inequalities) is equivalent to the 
fact that no point of the process N e defined by (|5.ip lies outside A = (0, z%\ X 
• • • x (0, Zn]. The latter probability equals exp{— Af tn ((0, oo) dn \ A)}, so 
that (|5.3p yields that 

P Mh) <zi,.. .,v(tn) < z n } = exp j-Emax (e**^/**) } ( 5 - 5 ) 

for all t\,...,t n e R and zx,...,z n € M. d . Applying this for n = 1, it 
is easily seen that condition (15.2p ensures that rj(t) is a.s. finite for every 
t £ R. Furthermore, the above argument shows that the finite-dimensional 
distributions of 77 uniquely determine the finite-dimensional distributions 
of N e . In particular, N e is stationary if and only if 77 is stationary. The 
following definition appears in [6j, however only for stochastic processes with 
values in the real line. 

Definition 5.1 (see |6j). A stochastic process £ satisfying f|5.2[) is called 
Brown-Resnick stationary if the process 77 defined by (|5.4p is stationary. 

Theorem 5.2. A stochastic process t € M, is Brown-Resnick station- 
ary if and only if 

¥>*!,...,*„ («i - id~ l l,u 2 , ...,u n ) = (pti+s,...,t n +s(ui ~ id~ 1 l,u 2 , ...,u n ) 

for all s,t\, . . . ,t n G R and u\,...,u n G R rf satisfying Y17j=i u l = ®> where 
1 is the vector with all components equal to one. 
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Proof. If we consider the additive system, i.e. let Xi = log(yi) and N(t) = 
{xi +£i(t),i > 1}, then the Brown-Resnick construction corresponds to the 
situation where {xj, i > 1} is a Poisson process on the line EI = {(x, . . . , x) : 
x E R} in R rf with intensity e~ x , x £ R. Then the result follows from 
Theorem ETT51 □ 

Since the measure A e is prescribed, the Brown-Resnick stationary pro- 
cesses form a subclass of stationary particle systems with special (namely 
diagonal) intensity measures supported by the line HI. Then £ H is the vector 
with all components being £ = d~ l X^f=i an d = £ ~~ £ H - The following 
result characterises all pairs of a Gaussian process and a Poisson process on 
H that yield stationary particle systems. 

Theorem 5.3. A Gaussian process £ ; such that (i>,£) is noi stationary 
for some v ^ H, and a locally finite measure A on i/te Zme EI generate a 
stationary particle system if and only if A is proportional to e^, 



PF(t) + 6(t) + c i/ A = , 
W(t) - \\a 2 (t) + c if A + 0, 



where W(t), t G R, is a centred univariate Gaussian process with stationary 
increments and variance o~ 2 (t), b is an additive univariate function, c 6 R 
is a constant, and £ IH ,£" L satisfy the conditions (ii)-(iv) of Theorem \4-l(>\ 

Proof. The sufficiency is easy to show. For the necessity, note that since A is 
concentrated on H, the projected particle system N v = {{v,Xi + > 1} 
is also a particle system generated by a non-stationary Gaussian process. 
The characterisation of univariate particle systems from [5] yields that the 
projected A is a exponential measure, and so is A itself, since it is sup- 
ported by one-dimensional subspace. The rest of the proof follows from 
Theorem HT6) □ 
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